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ABSTRACT 

The general objective of this project was to perform research In large 
adaptive antenna arrays for space technology applications. Specifically two 
tasks were considered. The first was a system design study for accurate de- 
termination of the positions and the frequencies of sources radiating from 
the earth's surface that could be used for the rapid location of people or 
vehicles In distress. This system design study led to a nonrlgld array about 
8 km In size with means for locating the array element positions, receiving 
signals from the earth and determining the source locations and frequencies 
of the transmitting sources. It was concluded that this system design Is 
feasible, and satisfies the desired objectives. 

The second task was an experiment to determine the largest earthbound 
array which could simulate a spaceborne experiment. It was determined that 
an 800 ft array would perform Indlstlngulshably in both locations and It Is 
estimated that one several times larger also would serve satisfactorily. In 
addition the power density spectrum of the phase difference fluctuations across 
a large array was measured. It was found that the spectrum falls off approxi- 
mately as f~5/2. 
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STUDY OF LARGE ADAPTIVE ARRAYS 
FOR SPACE TECHNOLOGY APPLICATIONS 


1. INTRODUCTION 

High angular measurenent accuracy and resolution implies large aper- 
ture. The resolution of an imaging device is in the order of the recipro- 
cal of the number of wavelengths across the aperture. Thus a five cm lens, 
which is 10^ optical wavelengths in diameter, has a potential resolving 
power of 10 ^ radians. There are many applications for which optical 
wavelength electromagnetic waves are not appropriate. Clouds and rain 

which can be penetrated easily by radio waves are essentially impenetrable 
to optical waves. 

Since microwaves are much longer than optical waves, a system which 
uses these waves to realize a resolution equal to that of an optical in- 
strument must possess an extremely large aperture. Most current radars 
operate with microwave wavelengths in the range from 3 cm (X-band) to 
30 cm (L-band). The size of a microwave aperture with a resolution of 
10 ^ radians would be 3 to 30 kms. 

Large aperture systems introduce fundamental problems. Microwave 
dishes larger than a few hundred feet are usually Impractical because of 
distortions resulting from gravitational and thermal stresses. Further- 
more, the silhouette of a large dish makes it impractical for certain en- 
vironments. The phased array is a legitimate approach to achieving a 
large aperture without imposing the mechanical and silhouette problems of 


the microwave dish. 
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Conventional phased array design, however. Is totally Inappropriate 
for a giant structure, for the following reasons: 

(1) Conventional phased array design requires that antenna elements 
be spaced on the order of one-half wavelength. Thus achievement 
of 10 ^ beamwldth would require more than 10^^ elements, which 
Is Intolerably costly. 

(11) Conventional phased array design requires a planar surface which 
Is Impractical for very large arrays. 

(Ill) When the aperture exceeds some size, which Is a function of the 
wavelength, the practicality of surveying the position of the 
parts of the array to the required 1/10 of a wavelength de- 

(Iv) creases rapidly. 

(Iv) Large structures, such as ships and aircraft and large apace 
vehicles are nonrlgld, requiring a dynamic measurement of and 
compensation for the time varying coordinates of each antenna 
element to an accurancy of 1/10 of a wavelength. 

A minimum approach to the solution of these fundamental problems Is to 
thin the array so that no more than the order of 1000 elements, say are 
required Independently of the size of the aperture. The elements are 
distributed randomly through the aperture so as to eliminate the grating 
lobes which normally would result from a thinned aperture. The array 
can be made to conform with whatever surface Is Involved. In the current 
design another order of magnitude of thinning can be accomplished because 
of certain restrictions In the demands upon the system. 

The work reported herein is a natural extension of the general re- 
search program of the Valley Forge Research Center. The central theme of 
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this work Is the development of techniques by which large, nonrigid antenna 
arrays may be used to operate in turbulent transmission media. The basic 
concept is that of reflecting radiated energy from a remote target, sens- 
ing the back-scattered radiation via an array of antenna elements which 
are scattered randomly over some surface or through some volume, then pro- 
cessing the received signals to form a high resolution image of the target. 
Such a system is referred to as a radio camera . 

The required technological innovations associated with the radio 
camera fall primarily in the area of adaptive operation of large, random, 
thin, conformal arrays. Fundamental research is being carried out in 
adaptive space-time signal processing; adaptive nulling; multipath supres- 
sion; optimization of image quality on the basis of feedback from the 

image; phase measurement techniques; angular measurement techniques; and 
signal parameter estimation in the presence of noise, clutter, and propa- 
gatl>m anomalies. 

Potential applications of the ladio camera are many: airspace, ocean 

surface, underwater or ground surveillance for military or civilian pur- 
poses (ground based, floating, airborne or satellite), high angular reso- 
lution imaging of internal organs in medicine, selsmlcs, and meteorological 
measurement . 

To this end, the NASA sponsored work (to be described in this report) 
was a first examination of the potential of the radio camera in space. 

In particular the design deals with a large spaceborne array with many 
elements randomly located in a near-planar surface for the formation of high 
angular resolution beams. The purposes of the array were: 
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(a) The location of a disabled air or surface vehicle in distress 
by direction finding on a microwave emergency emission from 
the vehicle. 

(b) Multiple frequency high angular resolution beam-forming for 
communicating with multiple ground sources and/or monitoring 
sources of radio interference. 

The work for this project was divided into two tasks as follows: 

Task 1: Design of self-cohering microwave system for high resolu- 

tion beamforming and location of signal sources. 

Task 2: Conduct a spatial correlation experiment to determine 

propagation limitations to the pi.vformance and useful size 
of a large, conformal array on the surface of the earth. 

Task 1 was concerned with identifying technology that may be useful in the 
application of large spaceborne as well as earthbound arrays. Task 2 was 
concerned with determination of spatial correlation of L-band signals re- 
ceived on the ground from a transmission from a geosynchronous satellite 
(ATS-6). 

For task 1, three system concepts have been developed. The first 
system concept was based on the use of fixed references on the earth's 
surface to synchronize a randomly distributed antenna array which is in 
geostationary orbit. Once synchronized, the array scans the region in the 
neighborhood of the synchronizing source to locate any < cher sources which 
may be present. Even though the concept proved to be theoretically and 
physically sound, its world-wide implementation was found to be impractical , 
since it required the ability to place beacons uniformly on the earth's 
surface which for geo-political reasons is unacceptable. Thus the second 


- 5 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 


system concept was developed. In thlr concept it is assumed that an 
element location subsystem can be created consisting of beacons clustered 
on the land masses and distributed over the earth. Given that assumption, 
it is only necessary that each array module in the spaceborne system be 
able to measure the phase of the radiation field at sufficiently frequent 
intervals and transmit this information to the earth. The complex samples 
and the locations from which the samples were taken constitute the only 
Information required by a signal processor to locate the angular coordinates 
of the source or sources. The "space" system essentially consists of three 
parts. The first is concerned with the surveying procedure to be used for 
locating the positions of the elements of the spaceborne array. The second 
part consists of the electrical hardware requirements, while the third part 
is concerned with actual location and identification of distress signals 
from the earth's surface. 

A surveying algorithm has been derived which is based on the use of 
phase measurements at the antenna elements of signals transmitted by a set 
of well determined randomly distributed beacons on the earth's surface. 

The element positions so computed are unique for the given set of phase 
measurements. The tolerances on element location, beacon positions and 
frequency have been calculated and found to completely satisfy the system 
specifications. It has been shown that the same algorithm may also be used 
for locating and identifying distress signals, thereby solving also the third 
part of the design problem. 

As for the electrical system, analysis was focused on the use of phase- 
locked loops at IF. The loop extracts the carrier from the modulated signal 
emitted by the caller and this is sent to each remote antenna which also 
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has Its own loop. The differential phase is measured at the Individual 
stations, quantized and sent back to the reference (base) station. The 
base station also extracts the modulation bearing the caller's address 
code. This Information Is sent via a downlink to the ground terminal where 
the data are recovered and processed by a resident computer. 

It has therefore been concluded that the second rystem design Is com- 
pletely feasible from the electrical system point of view. 

Finally a third system concept, called the Simultaneous Transform 
Array, was developed. The new array system operates In a manner analogous 
to an optical lens. The array elements phase-shift an Impinging target 
wavefront so as to focus it on an energy-detecting Image surface. The 
entire target image will be recovered from the image surface rather than 
through the point-by-point manipulation of array element outputs. This 
system was not developed analytically because It was beyond the scope of 
the present effort. Although the second system concept has been found 
satisfactory. It Is felt that this new concept should also be considered 
as an alternative design if needed. 

In Task 2, an experiment was conducted with the primary intention of 
obtaining some measured estimate of the maximum size of a ground based 
antenna array which could simulate a spaceborne system under different 
atmospheric conditions. The maximum baseline In the experiment was 800 
feet. It was found that the RMS fluctuations In the phase difference for 
all measurements (of signals transmitted by the ATS-6 satellite) as a 
function of spacing between receivers was no more than 1°. Hence It was 
concluded that an earthbound experiment with an aperture of 800 feet would 
be Indistinguishable from one Is space. 
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Thus In Section 2 the vsrlous system concepts sre developed and their 
respective performances analyzed. In Section 3 the results of the ex- 
periment are analyzed. Tn f action 4 It Is concluded that system concept 
two Is completely feasible from the electrical point of view, and It satisfies 
all of the system requirements. Further the experimental part of the project 
was satisfactorily completed . Finally, In Section 5 recommendations for 
further work are detailed. 

2. ARRAY- IN- SPACE SYSTEM CONCEPT DEVELOPMENT 

2.1 GENERAL OBJECTIVES 

The main object Is to design a large spaceborne array with many elements 
randomly located In a near-planar surface, for the formation of high angular 
resolution beams. The purposes of the array may Include: 

a. The location of a disabled air or surface vehicle In distress 
by direction finding on an emergency microwave emission from 
the vehicle. 

b. Multiple-frequency high angular resolution beam-forming for 
communicating with multiple ground sources and/or monitoring 
sources of radio Interference. 

In general we confine our attention to a spaceborne receiving array, 
thus avoiding the severe microwave power generation and distribution prob- 
lems that would result if the array were required to transmit rather than, 
or as well as, receive. While a geostationary orbit satellite Is concep- 
tually attractive, lower orbit cases may also be considered, since the 
extreme distance of the geostationary orbit may be prohibitive for some 


applications. 
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In line with the general objectives discussed above, it is possible 
to postulate a list of system characteristics that would appear to be 
necessary. 

1. To achieve beam formation from a sriceborn.i array giving resolu- 
tion of about 1 km on the ground would require an array size of 
H wavelengths where H is the satellite height in km. Thus at 
L-band, a geosynchronous satellite array would have to be 7 or 
8 km in diameter to provide the indicated resolution 6 (»60X/0 

degrees, when D is the array diameter). 

2. Radio Camera operation is possible since fixed or floating refer* 

ence sources can be provided on the earth's surface for synchro- 
nization of the spaccborne array as well as for calibration of its 
pointing directions. It should be easier to locate the position 
of an unknown earth-bound source with respect to a nearby (set 
of) reference source(s) than if absolute accuracy were required. 

3. Array elements can have" reasonable" gain since only the visible 
earth surface need be illuminated. Reference sources can be made 
strong enough (ground-based power and directivity) to insure high 
Si'IR at each element to permit accurate beam formation. Subsequent 
scanning of high resolution beams (in both space and frequency) 
then can permit adequate detection and hence location of unknown 
sources . 

4. Assuming 2000 elements located at random over an 8 km circle, the average 
spacing between elements will be on the order of 150-200 meters. 

Signal flow problems must be solved for RF, local oscillator, IF, 
video and power supply as well as position stabilization/calibration 
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of .leaenc structures. To collect dsts at a central point In the 

3 8 

array will require delays up to 4 x 10 /3 x 10 or 13.3 microseconds. 
The ability to achieve coherent spatially correlated beam formation 
would thus apparently impose the requirement that signal bandwidth 
be limited to no more than a few tens of kilohertz unless delay-lines 
were included in the design. 

5. It will be important to quantify how accurately it will be possible 
to determine element location. This will serve to limit the extent 
to which bean scanning can be performed and hence will determine 
required spacing of reference sources (in both position and frequency). 

6. L-band is probably a desirable design objective; lower frequencies 
would require greater size of array for comparable resolution while 
higher frequencies may make the system more vulnerable to weather 
and atmospheric anomalies although the smaller antenna sizes required 
make this an alternative worth considering. 

7. The geosynchronous orbit assumed in the above discussion provides 
continual satellite presence and simplifies satellite location from 
the ground. However, a lower satellite orbit can provide orders of 
magnitude more SNR for received signals as well as requiring much 
smaller antenna size for equivalent resolution/accuracy on the ground. 

It is apparent from the above list that while the concept shows promise, 
each item immediately suggests a host of questions that must be answered. 
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2.2 FIRST SYSTEM CONCEPT: Self-Cohering Redlo Cauere* 

The ala la to develop a ayatev for accurate dcterainatlon of the 

poaltiona and the frequenclea of the radiating aourcea which are lor.ired 

on (or cloae to) the earth's surface and whose radiation frequencies range 

between some minlmua frequency f . and soae aaxlmua frequency f . The 

Bln aax 

Idea Is to employ a large random receive array In a geostationary orbit; 
spectral analysis of the received signals and beam foraatlon after that will 
give the frequencies of the radiating sources as well as their angular 
locations within the array coordinate aystea; the longitudes and the lati- 
tudes of the radiating sources then can be easily found by relating the 
angular positions of the points on the earth's surface In the earth's co- 
or''l.-ate system to the angular locations of the source points In the array 
• .oordinate system. (See Appendix I) 

Our goal Is to achieve accuracy In position determination on the order 
of one kilometer. This can be facilitated if the beams on the earth's 
surface are on the order of 1 kilometer wide. Taking Into account that the 
distance between the geostationary orbit and the earth's surface Is roughly 
40,000 kilometers one finds that the required angular beamwldth would be 

Y - 25 X 10 ^ radians. (2.1) 


*This system, although technically feasible. Is discarded as Impractical 
because of the nature of the earth-based beacon subsystem which It requires. 
The reader may choose to go directly to the second system concept found In 
Section 2.3 which does not have the same limitation. 
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Thls Implies chst we have to deal with an array of size equal to 

D - - - 40,000 X (2.2) 

Y 

Where X la the wavelength of the received radiation. 

Although In principle we can scan throughout the whole visible (from 
the array) part of the earth's surface, scanning in the neighborhood of 
the edges of the earth's disc would not be practical because (a) the angle 
between the line of sight and the earth's surface becomes too small to 
permit satisfactory accurate position iQcation and (b) many of the radiating 
sources located in these regions will be screened out by mountains, forests 
and other obstacles. For these reasons, we will be concerned only with that 
part of the visible region for which the angle between the line of sight and 
the earth's surface is greater than 20°. The associated geometry is presented 
in Figure 2.1. 

As seen, the maximum scanning angle 9^ (measured from the line 
which connects the center of the earth with the location of the array) 

would equal 8.16°. Further, one can easily find (from “ 61.84°) 

that about 26Z of the earth's burface would be inside the region of inter- 
est . 

Determination of the positions of the microwave sources requires scan- 
ning in the (u-v)-plane on the earth's surface. Determination of the 
frequencies of the radiating sources will require scanning in frequency. 

If we generalize the meaning of the term "position" so that it includes a 

frequency coordinate besides the spatial coordinates, we can view the posi- 

tions of any radiating source as a point in the u-v-f space. Figure 2.2 


^^ALity 
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Figure 2.1 BASIC COVERAGE GEOMETRY 


shows the shspe of the volume within which the rsdlstlng sources of Inter- 
est will be located and through which we will have to scan. In the figure 

u and V are angular coordinates of the earth's center. In the special 
0 o 

case when the array Is horizontal, these coordinates equal zero. and 

f are the minimum and the maximum frequencies of Interest, respectively, 
max 
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The horizontal cross-sect ion of the 

for small values of u and v (such 

o o 

tem) and is an exact circle when u 

o 

equals 


volume presented is roughly circular 

as can be expected regarding our sys- 

■ V *0. The radius of this circle 
o 


r 


b 


sin6^ 


amax 


0.142 


(2.3) 



FIGURE 2.2 VOLUME OF SCANNING IN u-v-f SPACE 


At this point it is convenient to determine the number of angle- 
frequency resolution cells the system will have to deal with. Clearly 
this number will depend on the resolving power of the system in space and 
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frequency and on Che dimensions of Che volume in Figure 2.2. We have already 
made assumpclons regarding che size of Che angular resolucion calls and 
Che special dimensions of Che volume. From chem one can deCermlne Che re- 
quired number of angular resoluClon cells. Thus, since Che area of Che 
3 dB InCersecClon of a beam in Che u-v space approximacely equals *nd 

since Che area of Che horizoncal inCersecCion of che volume in Figure 2*2 
2 

is ve find che number of angular resoluCion cells Co be 

ft 

L - -2^ * 1.013 X 10® (2.4) 

e 


To find Che number of frequency resolucion cells we muse assume cercaln 
values for Che difference f - f . and Che required resolucion in 
frequency, 6. Arbicrarily we choose 320 MHz for (f^ ^^ - 
for 6 and for chese values we obCain Che number of frequency resolucion 
cells CO be 


Lj - - 3200 


(2.5) 


The values of (f - f . ) and 6 can be in principle larger or smaller, 
max min 

However chere is an upper bound on Che value of 6 above which we cannoC 
go. Namely, if we are Co achieve an angular accuracy approximacely equal 
Co Che widch of Che beam, we cannoC Cake for 5 a value greacer chan Che 
'array bandwidch. (This is discussed in Appendix II). 

Table 2.1 gives Che array bandwidch values Chac correspond Co a 40000 
wavelengCh array for four differenC working frequencies and aC Chree dif- 
ferenC scanning angles. (cbe angles are measured from Che broadside direc- 
Clon). (See Appendix IV). For che sysCem of concern Che values in che lefC 
column can be considered as an upper limiC of 6. 
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i 

1 

♦ 

^ 

frequency * array size 

, array bandwidth [kHz] | 

0 -2.5° 
s 

0 -5" 

s 

> 

0 

1 

• 

<x> 


GHz 

' 12 km 

489 

245 

i 1 

! 123 



L-band | 



! 1 

2 

1 

GHz 

1 6 km 

1 

978 

489 

246 I 

u 

GHz 

1 

' 3 km 

1956 

979 

491 



C-band | 




8 

GHz 

1 1 . 5 km 

3911 

195« 

, 983 

1 

...i. ■ — 


TABLE 2.1 ARRAY BANDWIDTH 


From (2.4) and (2.5) it follows that the total number of angle-frequency 
resolution cells is 

L - X Lj - 3.242 X 10^^ (2.5) 

To interrogate a volume in the u-v-f space we synchronize the array 

by a reference source positioned at some point (u , v , f ) , usually in 

o o o 

the central region of the volume of interest, and then scan around that 

point. How far one can go from the synchronizing point (u , v , f ), 

o o o 

or more precisely, what will be the shape and the size of the volume that 

can be covered by scanning, will depend on the. uncertainties of the array 

element location. We have shown in Appendix III, that in the special case 

when the rms tincertainties of the array element location are equal along 

all the three axes of a rectangular coordinate system, l.e., G - G ■ G , 

X y z 

the shifts in angle and frequency (Au, Av, and Af) from the reference source 
must satisfy the inequality 


o 


(Au)^ + 


(Av)^ 5. 


^lOf o ‘ 
o r 


(2.7) 
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In Che equation above f is Che frequency of Che reference source, C le 
the speed of light and - (o^ + i» the rms uncertainty of the 

array element location. It is easy Co see Chat the volume within which 
one can scan is an ellipsoid centered around the reference source posi- 
tion. The dimensions of this ellipsoid are given in Figure 2.3. As can be 
seen, they are inversely proportional to the rms uncertainty o^. In Che 
discussion to follow we will use the above mentioned special case as a 
model for analysis of oiir system. 

For an array as large as 40000 wavelengths it seems reasonable to 
assume an rms uncertainty whose magnitude is of the order of 10 wavelengths. 



FTCURE 2.3 MAXIMUM SIZE OF SCANNED VOLUME 
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If we make such an assumption, then the dimensions of the ellipsoid in 
Figure 2.3 will be such smaller than the dimensions of the voluise of interest 
in Figure 2.2. Therefore many reference sources will be needed to cover it. 
To determine the number of reference sources required, we will have to 
assume certain distribution of the reference sources throughout the volume 
of interest (Figure 2.2) and to find the required distances in frequency 
and space between the neighboring reference sources. For reasons of 
practicality we assume a distribution as the one given in Figure 2.4. Such 
a distribution can be obtai.ied if we distribute on the earth surface 
multiple frequency transmitters in such a way that in the u - v plane 
they form an approximately square lattice as shown in Figure 2.5. 


WOLIlPlt 

F.ltQUEKCy 
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FIGURE 2.4 SOURCE DISTRIBUTION IN u-v-f SPACE 
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FIGURE 2.5 SOURCE DISTRIBUTION IN u-v PLANE 


For the sake of simplicity in processing and to avoid unnecessary 
overlapping between the regions of two reference sources, instead of scan- 
ning throughout an ellipsoid we will associate with each reference source 
a scanning volume of the shape of a rectangular parallelepiped. This scan- 
ning volume will have orientation as that shown in Figure 2.6. Its dimensions 
2AUo» 2 Av^ and 2Af^ will be chosen in such a way that they satisfy 



■ Au 

o 


Av 

o 


( 2 . 8 ) 


o 




(av : 

' o'' 'lOf o ' 

o r 


(2.9) 


Note that this is the largest rectangular parallelepiped that can be 
circumscribed by the ellipsoid depicted by (2.7). 
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figure 2.6 ACTUAL SCANNED VOLUME 


Solving (2.S) and (2.9) one finds 


Af r 

o . . c 

- — ■ Au ■ Av Tn7~~r~ 

i o o lOf 0 

o or 


( 2 . 10 ) 


It is easy to see that in order to cover the volume of Interest 
(Fl{ture 2.2) by associating with each reference source (Figure 2.4) a scanning 
volume of the form, size and orientation given in Figure 2.6, the separa- 
tions between the neighboring sources (Figure 2.4) will have to be 


along the frequency axis: 

* 

in the u - V plane : 


a- 2Af 

o 

b- 2Au //*“ 
o 

- 2Av //2” 
o 


( 2 . 11 ) 

( 2 . 12 ) 


*We take b ■ 2Au / /T~ ■ 2Av / /T~ instead of taking b» 2Au ■ 2Av 
o o o o 

since we want to keep the orientation of the scanning volume fixed for 

the u - V plane, i.e., Independent of the orientation of the lattice of 

transmitters in the same plane. 




- 20 - 


ORKjINAL paqc is 
OF POOR QUALITY 




•/ 


Coahlnlng (2.11) und (2.12), And iiHing (2»10) yield* 
C 



(2.13) 


b 


C 

5£ o“/r 
o r 


(2.14) 


In Table 2.2 we give the values of a and b for several values of f . 

o 

assuming 

a - 20X - 20C/f 

r o o 


- 20X - 20c/ f ; R - 37000 km 
'TO o 



a 

T ■■ 

b 

' bxR 

2 

20 MHz 

• 

0.007 

i 261 km 

4 

40 MHz 

A 

0.007 

261 km 

8 

80 MHz 

1 

j A 

0.007 

261 km 


TABLE 2.2 MAXIMUM SEPARATIONS IN SPACE (b) 
AND IN FREQUENCY (a) BETWEEN THE 
REFERENCE SOURCES. 


From (2 >14) we calculate the number of transmitters required 


2 2 2 2 

irr‘ 50Trr^o^f‘ 

N k . 

T 2 2 

b C 


The number of discrete CW signals p<>r transmitter follows from 


M - 


^max ^mln _ 320x10^ 
a (C/5o^) 


(2.15) 

(2.13) 


a 


(2.16) 
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Tahla 2.3 gives Che values of and M for several values o£ f^, aasua- 
Ing - 201^ 


0 ■ 

r 

20X - 20C/f 

0 o 

1 

— t 

f^lGHzJ 

*^T ■ 

} 

M 

2 

1267 

16 

4 

1267 

8 

8 

1267 

4 


TABLE 2.3 NUMBER OF REFERENCE TRANSMITTERS (N ) 
AND CW SIGNALS PER TRANSMITTER (M) * 
REQUIRED. 


NOTE: 

Since the spatial separation between the reference sources is taken 
to be by /2 smaller than the spatial dimensions of the scanning volume 
(rectangular parallelepiped), there will be overlapping between the scan- 
ning regions associated with two neighboring reference transmitters. 
Therefore the actual number of angle- frequency resolution cells to be 
interrogated will be by a factor of 2 larger than the value given in (2.6). 
Thus we have 

L ~ 6.482 X 10^^ (angle-frequency resolution cells) 

In conclusion, the system analysed above has the following disadvantages: 
(1) The high beacon density Imposes a requirement that must be distri- 
buted uniformly anywhere on the earth's surface. This is Impractical since 
beacons cannot be located on the oceans and they can only be located in 
countries willing to cooperate. 


t # 
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(II) A large nu^er of coaputaClons Bust ba aadc for the large miobar 
of calla (In the u-v-f apare) that nuat be aearched. 

(III) The large nunber of beacona, and the requireaent that each beacon 
has to tranaalt at a nuaiber of frequenciea (Table 2.3), Impoaea a major 
Inatallatlon and maintenance problem. 

For theae reaaona It waa decided that thla ayatem concept ahould be 
abandoned for the more feaalble dealgn of ayatem concept 2, tdilch la described 
next. 

2.3 SECOND-SYSTEM CONCEPT: Self-Surveying Radio Camera 

The ayatem analyzed above Imposes an unacceptably small spacing 
between beacons. It Is not the large number of beacons per se that repre- 
sents the problem, for the amount of equipment required for constructing and 
maintaining these beacons Is miniscule compared to the spacebome part of 
the system. Instead the uniformly high beacon density imposes a requirement that 

beacons may be placed arbitrarily anywhere on the earth. Beacons can not be 
so deployed, however. They can not be conveniently located on the oceans and 
they can only be located in those countries willing cooperate. Thus the 
high beacon density is an intrinsic limitation. 

Another factor which makes the high beacon density system questionable 
is the suspicion that with so large a number of beacons on the earth it should 
be possible to dynamically survey the element positions of the spacebome 
array and thereby preclude the necessity for self cohering. 

With these factors in mind an alternate system concept was developed. 

The concept is based on the assumption that an element location substation can 
be created consisting of beacons clustered on the land masses and distributed 
over the small angle subtended at the geosynchronous satellite by the earth. 
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Glven chat assumption It is only necessary Chat each array module In the 
spaceborne system be able to measure the phase, or at most Che amplitude 
and phase, of the radiation field at sufficiently frequent time intervals 
and transmit this information to the earth. The complex samples and the 
locations from which the samples were taken, constitute the only Information 
required by a signal processor to locate the angular coordinateii of the 
source or sources. 

Thus in Section (2.3.1) the problem of surveying the elements of the 
array is discussed. In Section (2.3.2) a position finding algorithm, for 
locating the transmitting sources, is introduced, which has many conanon 
features with the surveying algorithm of Section (2.3.1). In Section (2.3.3) 
Che required electronic system is discussed. 

2.3.1 ARRAY ELEMENTS SURVEY 

As a consequence of the study of the sec'*’'J system concept, it was 
found imperative that the positions of the elements of the randomly dis- 
tributed spaceborne antenna array should be surveyed [1]. Thus given a 
set of transmitting beacons on the eart I's surface, the objective is to 
locate the position of the elements of the array using a set of phase 
measurements of the originals transmitted by the beacons together with the 
knowledge of the transmitting frequencies and the positions of the beacons. 
The algorithm of solution is based on the use of phase measurements at the 
antenna elements of signals transmitted by a set of well determined randomly 
distributed sources (beacons) on the earth's surface. The results so de- 
rived are unique for the given set of measurements. The algorithm is essen- 
tially based on the method developed in [21, (see Appendix V), for surveying 
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•n Array In the far field. In thla aectlon, we will sioq>ly aketch the 
outline of the algorithm and briefly atate the reaulta for the apaceborne 
array. A full deacription of the algorithm and the main reaulta may be 
found in Appendix V. 

Consider a set of N receiving elements which make up the apaceborne antenna 

T 

array. Let jr^ ■ (x^, k^) be the vector of unknotm rectangular coordi- 

nates of the ith element, with resoect to a reference eleisent at the origin. 
Assume that for sinusoidal excitations induced in the dements, the phase 
at each element relative to reference element at the origin can be measured. 
The objective is to compute jr^^. 

Suppose a beacon (on the earth's surface) , with known coordinates, 
radiates energy at wavelength X, arriving at the array of elements as a 
plane wave at elevation (90 - azimuth . We will characterize 

such a beacon by its wave vector 


kj « ^ (sin cos , sin sin ^ , cos 9^) 


0 . 17 ) 


“ X 

where a. “ sin 0 . cos <t>, b. ■ sin 0 . sin <{),, and c c cosO.. The absolute 
j J J J J J J 

phase at the ith-elensnt of the array, relative to a reference element, 

, due to a signal from is given by 


2n 


♦ij ■ X ‘“j*! 


(2.18) 


The sueasured phase is in fact * ('i'^^j) (rood 2 tt), i.e. 
^ij ■ X<“j*i * "j>'i ^ =j'i’ * 


(2.19) 


j - 1. 2, 


a • 


■A 
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1 


whnti* n^j 1 m huow InttfKMr, nnd M !• th« munhi>r of b#acon«. Bquactm (2.19) !• 
ri-pratod for rach clcnwot of thu array, 1 - 1, 2, .... N. For tha ith-alcmane. 
(2.19) may ba ramrlttan in matrix form aa 




il 


C- 


IM 




"2 S 


"n '’m ■'h 


L*l J 


11 


12 


IM 


( 2 . 20 ) 


l(*r M 2 . «iul aMiiutalny, flu* t)oacon» aro randomly dlstrlbuli‘d. It la aauy to abow 

that tha nulutlon vactor r^ ol (2.20) la aaaantlally unlqua [2]. Thia raault 
ni.iy tluMt ho uiu'd to di t orn’tno , 1 1 om t hi* vi'i'tor of olmarvod phaaoR , tha 1th- 
«*lpr.ii’nl posit Ion voctor 

In I'.rnoiol tlio ||•lallvl* piiato* voolor nii*aMuri*d at the array la “ 

- - ' V .j. * ‘ 

1'.',,. 0 .... .whara ^ In iji corrupt od hy random phase 

Jlia^ li'l I 1 

orroiH. A least Kt|uorru lit ij to tha observed phase vector Jij therefore 
piovides one estimate ol r^. Koi I'als purpose, an error function F(r^) In 
defined as follows; 


M , 2 

K(r ) ~ i. Uii’, ,(r.) • 4/ l(mod*?n)} 
. 1-1 ^ 


( 2 . 21 ) 


w’a *re (\>mod*?n ^ x 1 2nn lo: t;ome Intep.er n such that -lt<x4-2itn<x , with 

V (t ) helnp, the relative plnise shift (modiila 2ti) generated at the Ith elemant 
1 J I A 

hy the ith beacon. We are lookin(t for a position wetor r^ minimizing the 

error function; clearly. In the noise free case whan - i|; , wo have 

^ A — 

r 1 .. fj* - the true piisltlon vector. | 

Ciiaslder Inp. the overall characteristics of tills error function. Initially 

assumlnK the nolsa free situation, whan ooc In thn vicinity of r^*, the 

qMatlty under tha Munmatlon In (2.21) actn as the square of a random variable 

uniformly distributed over (-n.ti). Thus the average value F of the error 
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fund Ion In this case is given by 


F 



( 2 . 22 ) 


whev is close to r^* we can obtain a first order approximation to F(£j^) by 
it.si;’,; i Taylor series expansion around r^ ■ —1*’ yields the expected 

value E{p(Ar^)), given by 

ElF(Arj)) - + FAy^^ + GAz^^] (2.23) 


where F and G are given by 

F«y( cos coa^ ^ (2.24) 

G ■ -7 I 1 + cos a + cos^ a 1 (2.25) 

3 u u 


This analysis of the error function is for the specific case when the sources 

are uniformly distributed over an area whose boundary subtends an angle about 

the z-axis. For the array in synchronous orbit, the angle of view is only 17* 

and hence a - 8.3*. 
u 

Using the above function, the uncertainty in element location was found to 
be of the order of 1/3X in the x-y direction and O.IX in *'.he z direction, assuming 
phase errors of + 20 * [ 21 . 

Tha tolerance on beacon position may be computed as follows. Consider 
the tollowinij figure, which rcpre.«»ents the worst case, when the 1 th element 
and the refert'ncc element are perpendicular to the direction of the signal from 
the beacon. We will assume that errors in measurement of phase due to noise 
and uncertainty regarding beacon position are within the range of i 20 * (a fig- 
ure vhiclt Is easily attained with present day technology). From the figure 


on the next page it is easy to see tnat the phase error 64 / 
position of the beacon is given oy 


ij 


due to change in 



2 n(r^A) 


< tolerance on phase measurement ■> A ^ 

assuming no noise •• 2L 

9 ^ 
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Assuming ■ 1 km and X ■ 50 cm, then A 0.28 x 10 ^ rad. In other 
words If the antenna array Is In geostationary orbit, at 40,000 km from 
the earth's surface, assuming X ■ 50 cm, r^ - 1 km and phase tolerance of 
then the tolerance on beacon position l.l km which Is very easy to 
achieve. 

The tolerance on the frequency of the transmitted signal may also be 
computed [2]. Assuming a transmission frequency of 1 GHZ, the tolerance 
on the frequency Is ± 4 MHz. This Implies that the beacons do not have to 
have very tight specifications. 

The algorithm used for computing the position of the elements of the 
array Is based on using a number of beacons which are near the center of 
the field of view (l.e., small aperture) to obtain a coarse grid on which 
to search for the element position. This Is then followed by a finer grid 
formed by adding a number of beacons at the edges of the field of view 
(l.e., large aperture). For this purpose, a study was made to locate the 
set of beacons for different areas of search. Thus: 
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(i) To cover the Atlantic; beacons are placed in Western Europe and 
the East Coast of U.S. 

(it) To cover the South Pacific: beacons on the West Coast of U.S., 

Hawaii and Australia. 

(iil) to cover the Indian Ocean; Australia, Islands in the Indian 
Ocean and South Africa. 

2.3.2 SOURCE POSITION LOCATION ALGORITHM 

The objective of the proposed space system is to estimate the angular 
coordinates of a radiating source on earth, based only on phase measurements. 

In this section we examine two different approaches that have been proposed 
and suggests a hybrid that possesses the advantages of both approaches. The 
desired accuracy is 25 prad corresponding to an error on the earth of 1km. 

A. Interferometry 

This is a classical technique Chat has been used in radio astronomy. Its 
use for azimuth angle estimation has been studied in a VFRC report [3] from 
which some analysis is borrowed. The method basically uses a series of 
antenna baselines of different lengths, in which the final angular accuracy is 
determined by the largest baseline and a shorter baseline is used to resolve 
the angular ambiguity in the next larger one. Some details are given In (4). 

-3 

If the largest baseline is D,, we need 6D„ “10 D„ where 6D„ is ®51. 

This accuracy cay be achieved if the antenna baseline is almost rigid. 
With a moving baseline, however , such an accuracy is doubtful. Nevertheless it has 
the ->dvantage of computing esticates in a shorter time period than the second 
method to be described now. 

B. .‘"'-lar location by Least Square Fit Method, (LSE) 

Tnis approach was reported in [5] and uses phase measurements from a ran- 
dom acra/. it is a brute-force method that involves scanning 

over rh:; whole region of interest in the angular space to get cilniraura square 
error. .Although th.e earth is within a 17* cone, the azimuth i could be any- 
where Lthin (0, 2 tt) and a search is too time consuming for the desired angular 

accur . However, the tolerances on element positions are more relaxed than 
for a:i interferometer, tlic working of which as mentioned above hinges on the 
accuracy of the largest baseline. 


C. Prjposcil M**thod 
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Both Mthoda A and B uro attractive In complementary waya an • it aeesa logi- 
cal to consider a combination. We sliall use two short, almost perpendicular, 
r.eni-rifiid booms for interfcroraeter baselines as shown in Figure 2,7, The raat 
of the array system will be discussed after examining the interferoaater. The 
zeroth (reference) element is situated on the satellite itself. Hie objective 
is to use the interferometer to give quick, rough angle estimates without facing 
the tolerance problem alluded to earlier, then switch over to the LSE method which 
uses fbe rough estimates as initial values to achieve the final eatimatas that 
meet the accuracy requirements. 

In this method, we make the following first-order assumptions: 

1) the satellite coordinates relative to earth are known 

2) the Incoming wavefront is planar 

3) there Is no multipath Interference 

4) negligible propagation/refractive effects along its path 

5) the elements do not move significantly during the phase measurements. 



I 


/ 




L 



L, CiT 'lOO X *4 

a' 6000 X 




FIGURE 2.7 SATELLITE ANTENNA SYSTE71 
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Using the Interferometer 

It is shown In [4] that the coar&c accuracy of 3-5 mllliradlans can be a- 
chieved with 2 baselines per Interferometer with 2A, “lOA to 15X. In the 
f.pace environment and with baselines of these magnitudes, it la assumed that 
variations due to temperature, gravity, etc. arc negligible. The procedure 
for determining the source angles 0^, In the satellite coordinates in 
Figure 2.8 are ahovn In [4]. Thaae vaimms are used aa Initial guesses for the 
next step, as follows. 



PTniRF 7 R A-N'GLE RELATIONSHIP BE-n/EEN AMTEIE.’A ELEMENT 
AND SOURCE DIRECTION. 

Using the LSE Method 

Tnc final accuracy in angle estimation is determined by the extent L of 
the array. The accuracy can be viewed in the context of pointing errors in an 
•array since, in that situation, the phase shifters are adjusted to point a 
beam lit a certain direction while in the situation here the phases induced by 
a source in a certain direction are being measured. Thus wc assume that for 
both elevation 0 and azimuth 

“o' % ' i 


rad 


(2.26) 
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where is the standard deviation of phase measurcncnt error, X/L is the noa- 

i< 

Inal array bcamwidth and N is the number of elements in the array (tha error in 

azimuth could be an order of magnitude higher if the elevation angle 0^ Is 

close to V because of an additional term l/cos9_ for o., 

2 T9 

Ihe number of elements N is to be large enough to suppress the 'side-lobes' 
in the error function described in [3 ] and for more accuracy. ^ probaballstic 
approach is attempted in [4] to determine the number needed. We 
then need to determine hox^ the elements are to be distributed. This is 
doi'e as follows : 

Given that the original uncertainty from the Interferometer is 3 mr, let 
the next order of accuracy to be achieved be 0.32nr. From (ft), with " 15, 


3 X 10'^ /l5 


Al33 assume a, 


“O. ^ + o. 

«>.t «>,q 


2 +o2 


9 ,el 


= 0.5 (rad) , where 


a. ^ (thermal 
9,t 


noise) =* 0.18 rad with SNR “15db, 0 ^ ^ (quantization) “ 0.06 rad with 

5 bits representation; and 0 ^ (element position) ^ 0.45 rad , corresponding 

to ar error of about Jq . Then, 400X. The algorithm begins with the 

ir.terfir?r. 2 ter estimates 9 , ^ and searches over the region 9 + 0.003. + 

00 °o— 

0.003 rad in steps of 0.3 mr, requiring 400 steps to get estimates 9^^, 

let the final uncertainty required be 25prad and let N 2 •• 20 elements. 

Then 1, =" 5000X. This Iteration requires “ 576 steps. The two Iter- 

ations take a total of about 1000 steps for^evi?uating F(9,<5),the error func- 
tion ir. [5]. 

7;.e antenna element distribution is Illustrated in Figure 2.7. if possible 
it is '.unstrained to lie almost in a circle rather th.in a sphere with the 
earth direction at its broadside. 
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2.3.3 ELECTRONIC SYSTEM CO.NSIDERATIONS 

The survey and source estimation algorithms have been discussed in Sections 
2.3.1 and 2.3.2. As is shown later the computational requirements are excessive. 
The survey itself takes M seconds (assuming 1 sec per beacon frequency and M bea- 
cons) Just to gather all the phase information and the elements would drift 
slowly even within this period. 

The algorithm calls for an updated survey on detecting an emission on earth. 
The satellite then switches back to that particular caller frequency band to 
collect the necessary phase information to calculate the caller's angular posi- 
tion. This is Immediately followed by a second survey. It is hoped that the 
element locations can be approximated by an 'interpolation' of the two survey 
results (which themselves are not very accurate owinj to the element drifts). 

All the information is sent down to earth and is stored on arrival in the 
computer. Note that the reference coordinates and the transformation parameters 
from satellite to earth coordinates at the appropriate times have to be known. 
Admittedly this is a weak link in the design because not only must these be ac- 
curately known but those values at an earlier time are required owing to propa- 
gation delays both in the satellite electronics and the transmission path of 
40000 km. 

Given that the above-mentioned information can be obtained, all the data pro- 
cessing need not be done in real time. 

Sources 

The sources on the earth are of low power, with low antenna gain, operat- 
ing in the L-band. The power-gain product is of the order of 10 watts, at most 
100 watts. Low power is assumed. The transmitter will be a small equipment 
easily handled and not requiring high performance. Low autenna gain is assumed 
so that care need not be taken by the user in orienting the antenna. The rad- 
iating frequency is controlled by a low cost crystal driving a frequency multi- 
plier which results in one of N allowed frequencies occupying a band of width 
W Hz. Each transmitter will have a simple means to impart a simple address 
code or "call letters" of the user. One such possibility is a series of toggle 
switches or buttons. A 20-bit code would permit 10 addresses or separate 
users. Simple ON-OFF keying or binary phase shift keying would be adequate. 

Information provided by Dr. Sajjad Durrani of NASA indicates that 15 MHz 
at each end of the 20 cm wavelength satellite band is relatively free at the 
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present tloe for such use. In the design discussed below, therefore, 
the value taken for W Is 30 MHz. 

Functional Block Diagram 

Figure 2.9a 1^® ® block diagram of the system. The right-hand section Is the 

! 

reference element; the left-hand section is the array element. The signal 

from the source is radiated at frequency f . The front ends of both equip- 

8 

nents are the same. Each has an antenna of conventional design, ouch as 
a dish or a horn. The antenna diameter is 2 to 3 wavelengths such that 
its beanwidth covers the entire earth. The RP bandwidth is W. A mixer 
is driven by a scanning frequency synthesizer. The IF bandwidth B*W/K is the 

spacing between the N allowed frequencies within the band W. B is deter- 
mined by the long term stability of the transmitting sources. Assuming 
a low cost oscillator, a stability of 1 part in 10^ is readily achievable. 
Permitting some margin for error, and assuming a wavelength of 20 centi- 
meters, B is of the order of 10^ Hz. 

This value is taken to be the drift tolerance of the source trans- 
mitter. It is not the short term bandwidth of the transmission, however. 

_2 

Instead the sources are assumed to be stable for the order of 10 seconds 
corresponding to a line spectral-width of 100 Hz. This fact is utilized 
by following the IF amplifier in the reference module by a bank of narrow- 
band filters each of 100 Hz bandwidth. A bank of 1000 filters is re- j 

quired for this process, probably implemented digitally. Envelope detec- 
tion and low pass filtering for 1 second follow. Following integration, 
threshold and decision logic determine whether or not a signal is present, 
extracts its call letters or address, measures the frequency of the source, 
possibly alters the search pattern of the frequency synthesizer, and de- 
livers signal-frequency information to the frequency synthesizer so that 
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A.HRAY MODULE REFEilENCE MODULE (IN SPACE STATION) 

. s 



FIGURE 2.9a 


BLOCK DI.AGRA.MS OF REFERENCE .\ND ARRAY MODULES 
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Che latter may provide Che proper signals to the array element modules. 

In the lower part of the figure the received signal Is heterodyned to the 

center of the IF passband where it Is coherently quadrature demodulated. 
The signals of both channels are Integrated In low pass filters of 100 Hz 

bandwidth. Amplitude and pliase measurements can then be made. Only the 

latter Is shown. The phase measurement Is smoothed In a 1 second time 

constant low pass filter, after which the analog voltage is sampled and 

converted to digits. Five-bit precision is adequate. The element address 

is added to the five-bit word. The address requires log^N^ bits where 

N is the number of elements, 
e 


SHR Analysis 

The signal-to-nolse power ratio in the IF in the reference receiver 
and in each of the array element receivers is given by 


AitR^kTBF 


(36) * 


In (36) is the transmit power source on earth, Is the source antenna 

gain, Is the area of the receiving antenna in each spaceborne module, 

R is the (geosynchronous) range (40,000 km), k Is Boltzmann's constant, 

T Is Che effective temperature of Che receiving systems, B Is receiver 

bandwidth (10^ Hz) and F Is the noise figure. Since the receiving antennas 

are painted toward the earth, a nominal earth temperature of 290° Kelvin 

-21 

is assumed, making kl - 4 x 10 wates/Hz. At L-band a noise figure of 
a few d3 is available. To provide some margin F 10 dB is assumed in 
this analysis. Assuming PG^ •» lOW, Aj^ “ 4X and X = 20 cm, (36) evaluates 
to -17 dB. Coherent integration in the filter bank of the reference re- 
ceiver and in the quadrature low pai s filters of Che array element re- 
ceivers narrows the noise bandwidth to 100 Hz, thereby increasing Che 


I 

J 
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slgr.al-to-noise ratio by 30 dB to 13 dB. Noncoherent sooothing by a fac- 
tor of 100 to 1 further increases the signal-to-noise ratio by nearly 
20 d3. (The increase in signal-to-noise ratio through the post detection 
filter is approximately linear with the amount of integration because 
the input signal-to-noise ratio to the post detection filter is so much 
greater than unity.) The result is in excess of 30 dB. 

The angle estimation is a function of the number of elements in the 
array which in turn will be Influenced by the angle-measuring or direction- 
finding algorithm eventually chosen for the system; in general the larger 
the number of elements in the array the smaller is the SNR requirement 
per element. A conservative bound is obtained by considering the two ele- 
ment interferometer. The available angular accuracy is given by [26] 

A0 “ -V (37) 

Lr^ 

In (37) L is the size of the interferometer base. The desired angular ac- 
curacy is the order of 1 km at a distance of 40,000 km, or 2.5 x 10~^ rad. 
Taking Che signal-to-noise ratio to be 30 dB the maximum required array 

size is L ^ ■ 1440 m. A still more conservative estimate is ob- 

max rmao 

tainad by assuming 30dB SNR as a margin for system and component troubles. 
For example a user in distress on the ocean cannot be expected to align 
the antenna correctly or even, perhaps, to ensure that it is not partially 
submerged. A 30dB margin would materially enhance that probability that 
the system works when it is needed. In this case the required array size 
becomes L = A/40 ■ 8 km. 

Both numbers are large. An array compatible with this system concept 
will be measured in kilometers. 
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Satellite Electronics System 

It Is proposed that the earth-borne transmitters emit bursts of 
CW and BPSK signals alternatively, the latter containing a header for 
synchronization followed by the code letters. The CW Is used for locking 
the loops as discussed shortly. BPSK Is also used In the downlink to 
the earth station where data processing is to be done. 

A more detailed system design Is shown In Figure 2.9. The front 
end Is modified as shown. The IF signal Is fed Into a Costas loop which 
conveniently provides both an unmodulated reference carrier r(t) for 
phase comparison as well as the demodulated binary code z(t) that Identi- 
fies the caller. The loop must have a bandwidth large enough to track 
the center frequency of the BPSK Input but small enough to avoid track- 
ing the modulation Itself. If Che loop filter Is chosen to have a cutoff 
at 50 Hz to acconmodate the 100 Hz oscillator coherence bandwidth expected, 


Chen we require 



TO MUX 


FIGURE 2.9 SATELLITE ELECTRONICS SYSTEM 
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the kuylng tu exceed 50 Hji., *'ay a 50vi hpa rate. (A conflict with the detcctioa 
circ.Ltry would arise if the input ia continuout BP3K becauaa the power is 
aprejvl over several detectors in the bank. Tliis ia the reason for introducinf 
shore bursts of CM between the modulated wavefom). The t'vo LPF's in the 
C-ost^s loop I-Q arms nust be wide enough to pass the BPSK modulation, say 1.2 
kHz. It is shorn in the appendix in (61 that v require 500 bps rate BPSK 
transmitters operating at about 4 watts output. 

• 

isBPF-envelope detector bank has JOOO channels and is designed to indi- 
cate i:t voich bund the center frequency lies. Assuming that the modulation is 
PN-lik-i, then the PSK bpectrum Is of the sinx/x form with the main lobe within 
^ li ^ Since the NBPF's have a bandwidth of lOOHz only, there would be 

little vnergy at their outputs even with u 1 sec signal duration. However, 
the o ::. ..t of one of the NBPF's would be large when the input is a CW burnt and 
this .• .:'d regi '.Lor a hit at the detection circuitry with the appropriate thres- 
hold. 

On rcco#;nicLnc a hit in one of the 1000 channels, the logic stores this 
infon.ition as well as the PvF band tj. anl switches over to the survey frequen- 
cy band for an update on the clement locations and sends the phase and clement 
addro.ss InfornuLlon via the data link. Tlie only difference between surveying 
and listening to a caller ia that the former is CW and no Idcntif icatlon code 
Is nv4*dt-l because each beacon has its own nsolBned frequency. The different 
procedur*:a in handling both situations can he incorporated in the logic to save 
time, rather than including it in software. 

On comjilet ion of the survey, the logic f.wltchos tlic RF local frequency back 
to fj^ and sends a control signal that adjusts the free running frequency of the. 
digital VCO in the Costas loop so tb.it the IF signal at would be acquired. 

Onci in lock, the loop provides the dc.iodulated binary bit stream rvt) for the 
soiirc" address and a coherent IF carrier r(t) to be piped to the remote elc- 
re-..rs for phase co.mpar Ison , bearing in mind that delays to every clement must 
bo equalized. Note that a filter matched to a bit in the code sequence z(t) is 
used to i.iaximizc the SNR and hence reduce the bit error probability. 

The other elements arc separated from the reference station and, in addi- 
tlo.i to the IF signal, mist he provided with (a) the RF local reference at fpg 
and fb) the digital control signal from the base station logic that allows thm 
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rcroce Costas loops to acquire the lnco:aIng modulated IF signal. This schera 
la herter than bringing the various IF signals back to the base station ba- 
cau le those signals are weaker and would cause degradation In tha phase com- 
parisons. On the other hand, an amplified reference r(t) can be sent out to 
st.ncions. However, this present scheme requires a ph.ase comparator In each re- 
mote terminal. Tlic results expressed In 5 bit words arc sent back together 
with the element address. The latter would not be needed If wire coimcunlcations 
is iis<. i. 

rack at the base reference, v.arlous Information as shown In Figure 2.9 Is 
asse.r.iicJ aa a packet preceded and followed by a preamble ond postamble re- 
spect ialy. This is sho»/n in Figure 2.10. n>c frequency Information (to the 
nearest lOOliz) Is needed in both the surveying and source location computations. 
Since is one of the 300 bands and lies in one of the 1000 filter band- 

width.^, at least l‘l bits of data are needed. No word synchronizer Is built In- 
to the ..atellilc because this inform.it Ion is of secondary importance and can 
he extracted after the source coordinates are cumputed. Moreover it saves on 
hardvaro on the satellite. Consequently most of the binary Information from 
the earth-borne source in Figure 2. 1C, including its header, is relayed to the 
ground stat ion and this exceeds the source address which is only 20 bits wide, 
.'•.ucli more allowance r.iist be provided since the position of the he.ider is unknown. 
About 50-f)0 bits, which covers two full lengths of header and source code, should 
ttufilce. bach element address can he coded Into ^ bits for N about 20. 

So the address and phase information requires 200 bits in nil, making a total 
of about 300 bits including 20 bits preamble and 10 bits postnmble of all zeros. 
It should be noted that many tradeoffs could be conducted to optimize the system. 
For example the choice of 1000 narrow band filters and envelope detectors may be 
Improved by using lO-coarse and 100-flne detectors. However, this was beyond 
th« scope of the present effort. 
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The fir it half can be alternate +l's fur bit s>tic and the second half is a 
barker sequence for vord s>Tic. 

The bit streams are clocked into the multiplexer at their own rates but the 
output is a single serial bit stream at rate f^. As in Figure 2. 11, this is low- 
pass filtered to reduce the required bandwidth and then biphase modulated with 
a subcarrier prior to dnvn-convcrsion at about 1.5GHz. 
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FIGURE 2.11 lX)W;aiNK SYSTDI 

Earth Station 

The requirements on the amounc of information to be transmitted to ground 
is minimal but it is crucial that tine delays be minimized. For this reason, 
highest priority must be assigned to this if any form of resource sharing Is 



FIGURE 2.12 EARTH STATION ELECTRONICS. 


required. Alternatively it could be a dedicated channel (like SCPC) which Is 
activated on demand. Ue shall leave the up/down link design aside since it de- 
pends 0.1 current technology. 

On do'.m-conversion, the bi-phase modulated signal is demodulated by a 
Costas loop which does both carrier recovery and demodulation. The baseband 
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digital data are bit synchronized with an early-late gate type clffemlt, ***»pl« 

and clocked into buffers and hence into the computer after the preamble has been 

recognized, as shown in Figure 2.12. The data bits ought to be transferred as soon as 

possible and the buffer Is primarily a temporary store in case of delays. A postamble 

stream of zeros signifies end of transmission. In the case of surveying the 
* / 

processor has to wait until all M sets of transmission have been received be- 
fore any calculations can be made. The computer can distinguish beacons from 
sources by their frequency fp. Calculations of the bit error probability of 
the do'.mllnk data transmission are given In [6]. 

Computational Requirements 

We shall assume that a large computer shared with other users is avail- 
able. But highest priority must be given to our application. A minicom- 
puter Is not as efficient because the accuracy requirements demand at least a 
20 bit representation of the angle information (about 10 ^ rad) while minicomputers 
are by definition at most 16 bits wide. Of course, double precision could be 
used on a mini but it would be too time consuming. Moreover, cycle times in 
a mini are usually larger than for a regular large computer. 

The computer is used to? 

(1) Survey and track the element positions as they drift in space 
due to gravitational forces, solar wind pressure, etc. This Involves (a) an 
initial search over the whole region after deployment, (b) periodic surveying 
to trarV. the elements, (c) locating elements just before computing the source 
location when a caller is detected based on pre-and post-surveys. The latter 
ensures up-to-date information of the elements. 

(2) Use the phase information to locate earth-bome sources in- 
volving (a) a search for azimuth and elevation in satellite coordinates, (b) 
conversion to earth .'oordinates. 

The surveying procedure is controlled by logic on board the satellite and 
is sumrjrized in the state diagram in Figure 2.13. cotrputer waits for 

appropriate information from the satellite for data processing. Its operation 
is pictured in Figure 2.13. 
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FIGURE DIAGRAM AT FIGURE 2.13(b)COMPUTER USAGE 

SATELLITE 

i 

The survey alt^orlthiD Is described in detail in [2], while the algorithm I 

for computing the azimuth and elevation of the caller is described in Section j 

(2.3.2). Figures 2.14 a and b suianarize both algorithms and are presented for 

purposes of estimating computation time. 

Since the time taken to transform 0,^ from satellite to earth coordinates would 

be much less than that required to search for 0 and <(), they are neglected here. 

For the implementation of the survey, we will consider the steady-state 
tracking mode, since the initial search ' k only made once (unless the element 
drift is severe). At first a search over +5X about the previous coordinate 
values for each of x,y and z is made with a grid size AXj^, Ayj^, and AZj^ of 2X. 

This grid implies that Fj^(s) has to be evaluated 125 times. Subsequently a 
grid of size y , j and Jq 'reapectively is taken. All of the above would 
amount to approximately 320 evaluations of Fj^(s), and based on previous estimates 
this would require about 40 seconds of computing time. 

The source location technique, using the procedure outlined in [2], takes 
about 1000[600 x 20] psec » 12 sec for N - 20 elements. 
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FIGURE 2.14(b) SOURCE ANGLE CALCULATION FLOW CHART 
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It is seen that these computation times are large* Moreover if in the 
survey each frequency step is 1 second In duration, there is a delay of M - 
15 seconds for all the phase information to reach the computer. Hence a sur- 
vey takes at least 55 seconds. The 15 seconds could be reduced if less inte- 
gration time is used in the phase measurenents on the satellite. A faster 
computer could be used. It is because of this problem that the scheme shovm 
in Figure 2.13 is used. 

Unresolved Problems 

1) the large amount of computation time required especially in surveying idiich 
exceeds delay and data transfer times. 

2) a''ailability of accurate earth-satellite coordinate relationship at the 
instants of time that the corresponding phase measurements are made. 

3) motion of the elements during the phase measurements. 

4) maintaining the orientation of the antennas so that they face the earth at 
all times. 

5) effect of multipath (which might degrade the system). 

6) effect of scintillation and refraction in the atmosphere. 

Conclusion 

The electronic and signal processing system is feasible. The more expensive 
items in the satellite are the frequency synthesizer, the phase comparators and the 
Costas loops which are assumed available at the respective operating frequencies 
at nominal prices. Power amplifiers are needed to pipe LO signals to the remote 
elements. 

The demand on bandwidth is minimal and the earth station signal processing 
hardware is of the same form as that in the satellite. 
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2.4 THE THIRD SYSTEM CONCEPT: THE SIMULTANEOUS TRANSFORM ARRAY 

In thl* aectlon a new array concept called the Sloultaneoua Tranafora Array 
la Introduced, which may be referred to aa the third ayatem concept for thla 
array in apace. Although the aecond ayatem concept haa been found aatlafactory, 
it ia felt that thia new concept ahould alao be conaidered aa a viable alternative 
which may lead to a better deaign. 

The new array system operates in a manner analogous to an optical 
lens or parabolic mirror [7]; cither of these devices will form an image 
of all the points in a region about their major axis of symmetry 
simultaneously rather than poJnt-by-point as in the case of beamforming 
and scanning arrays. The optical lens or mirror introduces the 
appropriate set of phase shifts (delays) to convert the planar wave- 
front from a far-field point source into the quadratically shaped 
vavefront required to image the point source at an arbitrary focal 
point. When this is accomplished, points in the vicinity of the 
distant source can be imaged on a surface placed in the vicinity of 
the focus. The simultaneous transform array will function in a similar 
fashion. Tlie array elements will phase-shift an impinging target 
wavefront so as to focus it on an energy-detecting image surface. The 
entire target image will be recovered from the image surface rather 
than through the point-by-point manipulation of array element outputs. 

Fleure 2.15 shows a schematic representation of how a microwave 
simultaneous imaging system might be implemented. Each of the array 
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FIGURE 2.15 SIMULTANEOUS TRANSFORM ARRAY 
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clemenCs can detect energy arriving from the target region. Ihe Indivi- 
dual array elements amplify the detected signals, phase shift then and 
then retransmit the signals toward the image surface shown at the left 
in Figure 2.15. If Che appropriate phase shifts have been introduced to 
focus the target regions on the image surface, the energy distribution 
across the surface will be a microwave image of the target region. 

The image surface might be made up of any of the devices which are 
suitable for detecting microwave energy such as gas-tube detectors, 
small dipoles, or fluorescent chemicals. (Much important work has been 
done at the Moore S;hool with gas-tube detectors ior microwave and 
millimeter wave hol'^graphy [8]. A major factor in the choice of an energy 
detecting device will be the manner in which the image is to be recovered, 
e.g., whether it is to be viewed optically or converted into electrical 
signals for subsequent pr»'<cessing. 

The image plane will contain a receiver at its center which will 
transmit to a central processor the signal from a point-source beacon 
in the target region. Tlte central processor will adjust each of the 
array element phase shifters to maximize the signal observed at that 
receiving point. Such a technique is a modification of the beamfoming 
procedures which have been studied at VFRC. 

Once the adaptive focusing is complete the simultaneous transform 
array can function in a manner analogous to that of a lens. An off- 
axis target will produce a v.’avcfront at the array that is slightly 
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tilted relative to that of the beacon. The tllced phase front will be 
transferred through the array modules and will produce a similar tilt 
in the wavefront transmitted to the image surface. Thus the target 
point will be displaced in the image surface from the reference beacon 
by the angular distance between the target and the adaptive focusing 
beacon in the target plane. All other targets will be similarly imaged. 

A remarkable feature of the simultaneous transform array is the 
freedom of positions of the array elements and the lack of required 
knowledge of these positions. A second advantage is the eliminaticr of Che 
need for an image-forming signal processor. The primary disadvantage 
is the need for an image surface. The ultimate merit of the simult- 
aneous array concept will depend upon the particular application. In 
the following section two space applications are considered. These 
applications are suggested as vehicles for exploration of the first 
order properties of the Simultaneous Transform Array (STA) concept. 

Both applications to be described employ an array in space. Ihe 
first is a somewhat "conventional" surveillance or measurement applica- 
tion; a large array is used to obtain a high resolution image of the 
earth for such purposes as resource studies, or volume cells in the 
atmosphere for meteorological purposes. The STA concept is independent 
of wavelen^h ; it would permit high resolution imagery using relatively 
long wavelengths such as L-band, or it could work at much higher fre- 
quencies according to the detailed application, since the positional 
distribution of the array elements would not be a performance factor. 

In this type of application, the array will perform the function of a 
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microwave mirror end Che image surface will be placed between the 
array and the target region. 

The second example treats the array as a microwave lens. In this 
form the STA would be useful as a deep space radio telescope. This 
example will introduce some other attractive characteristics of the STA 
conct It . 


2.4.1 EARTH-VIEWING RADAR 

Consider an array of elements distributed randomly over a region of 
space; a lenticular cloud of array elements situated In a synchronous 
orbit 40,000 km from the earth, as illustrated in Figure 2.16 is an 
example. Assume that a beacon is placed on Che surface of the earth 
to serve as a beanforming point. The signals received at the array 
elements are cophased, forming a radiation pattern having a main lobe 
ir. t’r.e direction of Che beamforning beacon. The cophasing is accomp- 
lished with the adaptive beamforning techniques developed at VFRC, 
which involve measurement of the phase angles at the array elements 
relative to a local phase reference when the array is illuminated by 
the earth beacon. 

(In Figure 2.16 , a planar phase front (A) is d^own at the array due 
to the beacon, implying that the Earth is in the far field of the array, 
which is not a necessary condition.) 

A second beamforming process is accomplished by the same array 
modules. A tr.nnsnit beam pattern is formed directed at a nearby 
ir.age surface. The surface contains an aiming point receiver. The 
receiver and the image surface might be located on a satellite situated 
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bptwven the array and the Karth. The S/N ratio at the almlnit-polnt 
receiver can be m.tde hl^h because the distance la short (perhaps 1 
to 10 km). Also the position of the receiver with respect to the array 
Is relatively stable. With these characteristics a transmit beasi- 
fonnlng phasefront (B) can be obtained by feeding buck from the almlng- 
pelnt receiver to the array modules phase correction Information to 
each array element; this procedure might require considerable time 
but since It Is required only occasionally , this does not appear to 
present a problen. 

Fach array element thus stores two phases, and ll't* 

latter being the phase correction which must be added to the local 
array-phase-reference to produce the quadratic phase front (B) required 
to focus the array op the near-fleld aiming point receiver. 

If the beattf orr.ltig beacon is now turned off and the Kurth Is 

1 1 lutninated bv a radar transmitter ns indicated in Figure 2,16, target 

ret loot Ions will he obtained. A strongly reflecting point on the Karth 

would produce the tilted phase front (C) shown in Figure 2.16. The 

array modules would measure phase angles Tlr** phase difference 

between the stored beanforming angle and the measured target angle 

; is A . A Will be computed .ind subtr.tcted i rom the transmit beam- 

lornlng angle ■, to form tire new target transmit phase. Tiro result 
1 H 

is a tilted quadratic transmit phase front (0) causing the transmit 
bean to image the target at a polrrt displaced f rr>m the aiming-point 
receiver. 

To observe and recover the target Ititage it is necessary to surround 
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the alialng-polnt receiver by .1 detection surface consisting of an array 
of energy detectors which tntiy be scanned or read out sequentially. 

The inuige surface might consist of an array of dipoles and detectors, 
gas tube RF detectors, or any other transducer which produces an output 
proportional to Impinging RF radiation. It should be noted that these 
detectors do not have to be particularly sensitive as they are not 
responding to the low level return from the Earth target; Instead they 
are being illuminated by the relatively near-by array elements which 
detect and amplify he low-level Earth signals to produce arbitrary 
power levels at the image surface. In this way the random array 
functions as a microwave mirror, slmultaxieously transforming the 
Earth-region into an image on the nearby image surface without the 
necessity of polnt-by-point scanning of the array beam pattern. 

(Geometric distortions will arise as in any other diffraction 
system; as the angular extent of the target region increases and/or 
the distance from the array to the image plane is reduced, the distortion 
will Increase. Increasing the distance between the image surface <ind 
the array decreases the distortion but increases the size of the image 
surface. Reducing the area of view or using several Earth beacons 
which could be selectively activated also reduces the distortion and 
permits a smaller image surface without loss of resolution. Thus th«?re 
is a number of system parameters which interrelate and which might bo 
used to reduce the distortion and/or limit the size of the image piano. 

Another solution might be to structure the array ns a collection of 
subarrays each of which is limited in size so that botli the target 
region and the inuigc surface are in Its f.ir field and the tianerait and 
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receive wavefronts are both linear. The quadratic transmit wavefront 
required of a high magnification array would be approximated by a 
scries of planar wavefronts, one from each subamy. In still another 
form a subsystem would measure the element locations, the knowledge 
of which would permit correction for geometric distortions. 

Clearly the STA concept Is very flexible. A study Is required 
to determine the best ways of utilizing the concept for particular 
applications. 

2.4.2 RADIO TELESCOPE FOR DEEP SPACE 

The simultaneous transform array concept also appears suited to 
the Implementation of a large space-bome radio telescope. A schematic 

representation of such a system is shown in Figure 2.17. The array is 
used as a lens to Image a region of small angular extent onto an Image 
surface located on the opposite side of the array. The image surface 
is assumed to be relatively small; consequently the angle of view of the 
array Is small. The image surface might be mounted on a space vehicle 
so chat the direction of observation of the array could be changed by 
movement of the image surface. The distance between the array and the 
iTJige surface might also be varied to control the magnification of the 
array. Geometric distortion will be much less significant in this 
application because the angle of view of the array will be snail. 

The array will be focused in a particular direction by using a movable 
beacon on the image-surface side of the array. It is proposed to 
accomplish focusing by using a beacon frequency that is much lower than 
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that used for imaging. Use of a low frequency will permit focusing the 

array at infinity even though the focusing beacon is relatively close 

to the array. The far field focusing boundary for the array Is 

2 

approximately equal to D /X where D is the size of the array and X is 
the operating wavelength [9]. llius by using a low beacon frequency for 
adaptive focusing, the far field boundary is effectively close to the 
array and a nearby beacon can be used for far field focusing. The 
use of different frequencies for reception, transmission, and focusing 
is yet another operating possibility which might be exploited in 
some applications. 

The viewing direction of the radio telescope will be around the axis 
from the adaptive focusing beacon to the array. A logical and practical 
location for the beacon is on the image surface. The telescope will be 
aimed by physically moving the beacon and the Image surface. The 
magnification of the array can be varied by adjusting the distance 
between the image surface and the array. Following selection of 
direction and magnification the array will then be focused by trans- 
mitting a low frequency signal from the beacon. Tills source will 
generate a wavefront which appears planar at the array. The array will 
use this wavefront as If It were arriving from a point source at 
infinity on the object side of the array. After the transmit wavefront 
from the array has been focused on the Image surface, the focusing 
beacon can be turned off and the array can be used to produce high 
resolution images of targets in the viewing direction. 
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In conclusion It should be noted thst the sloultaneous Trsnsforn Array is 
different from any of the array systems %/hich have been studied at Valley 
Forge Research Center, and it carries with it a new set of problems. As the 
examples given have suggested it has a number of important features which 
warrant examination. One of the most important is the concept of high resolu- 
tion microwave imaging without knowledge of the locations of the array 
elements. This possibility alone would make the STA an important study area. 
Geometric distortion, measurement tolerances and system organizations need to 
be examined theoretically and tested experimentally. 
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3. ATS-6 TO EARTH SPATIAL CORRELATION EXPERIMENT AT L-BAND 

3.1 THE EXPERIMENT 

The prinary intent of the experliaent was to obtain some measured estimate 
of the maximum size of a ground-based antenna array which could simulate space- 
borne equipment. An array in space can be assumed to be in a transparent and 
stable propagation medium. Signals from the earth after passage through the 
atmosphere and ionosphere travel undistorted in a free-space environment to a 
geosynchronous array. On the other hand, the atmosphere and ionosphere are on 
the receiving end of a circuit from geosynchronous orbit to earth. Therefore 
the medium effects upon signals received on a large downlink aperture would be 
more severe than if the large array were in space. Experimental confirmation 
of any space system design must be performed on earth, however. Hence it is 
necessary to have some idea of the largest ground-based array which will pro- 
perly simulate a spaceborne system, under different atmospheric conditions. 

Since operation of the antenna array is critically dependent upon the 
spatial coherence of the phase of the revived wave across the aperture, the 
experiment was set up to gather some information about this function. An 
equivalent quantity is the variance of the phase difference as a function of 
spacing between receivers. It was this quantity which was measured. The maxi- 
mum baseline available for the experiment was 800 feet. It was found that the 
RMS fluctuation in the phase difference for all measurements was no more than 
the order of 1°, which is 1 to l^j orders of magnitude below tolerance. Hence 
it was concluded that an earth bound experiment with an aperture of 800 feet 
would be indistinguishable from one in space. It was also concluded that an 
array many times larger also could properly simulate a spaceborne experiment. 

Tne second measurement was the power density spectrum of the differential 
phase fluctuations. Although this function bears no relation to the rationale 
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for conducting the experlnent, it was possible to gather such information and 

to make the necessary calculations. It was found that the power density spectrum 

-5/2 

falls off with frequency approximately as f 

The experiment was performed by measuring the phase difference between 
two receivers tuned to the ATS-6 geostationary satellite L-band CW emission. 
Specifically, the objective was to determine the differential phase as a 
function of distance separation between the receivers. A strongly turbulent 
mediua would introduce large scintillations on the wave and the effect would 
be seen in the results as a large variance in the differential phase. 

The VFRC upper site was chosen because the satellite path would be in line- 
of-sigl*.t of the receivers over the entire period when the satellite was in 
trar.J^i: from approximately S70'*r to its final position of S 72*W, as shown 
in Figure 3 . 1 . 

Tr.il receiver setup used is shown in Figure 3. 2. Specifications are listed 
belo”: 

1. Low noise amplifier (LNA) lOMHz BW 

2. Mixer from TACAI! gear: (a) 3-section filter on LO input tuned to 

1503 14MKz BW; (b) 2-scction filter at RF input tuned to 

1550 fab:, 20 MHz BW: (c) A2 MHz output, A tiHz BW. Tnere is an 
IF preamplifier stage also. 

3. IF amplifier from TACAIi gear: all tube, A2 MHz, AMHz bandwidth. 

A. In-house line driver: isolator and amplifier for 50n cable driving 
for the IF output. 

5. Local oscillator FXR Inc. Model L772A: Measured noise figure, 

F = A dB; overall system BW = A MHz. 

6. A ft. dish antenna, right circular polarized, from MARAD. 
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PIGURE 3,2 RECEIVER SYSTEM 
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Flgur* 3«3 shows the phase measureaenC setup used. For separations 
of 40 feet or less, no power amplifier la required for the LO signal* 

For larger baselines it is essential because of cable losses at 1508 MHz 



FIGURE 3.3 EXPERIMENTAL SETUP 


The cables transmitting back the 42 MHz signals for phase comparisons were 
standard coaxial types (about 1.4 db loss per 100 ft). Special low loss cables 
(5 db/100 ft) were used to pipe the LO even with the power amplifier present. 
Owing to the absence of a power splitter at the time of the experiment, only a 
simple coax tee was used (which unavoidably caused some power loss). The phase 
meter was a HP 8405A vector voltmeter with a IKHz bandwidth and an accuracy of 
about 0.1*. The result was recorded on a paper strip recorder which had about 
0.3 Hz bandwidth. 

The experiment was conducted mainly at ground level at VFRC for the large 
separations. The roof top gave a better view of the satellite, but was useful 
only for separations of less than 100 ft. The actual receiver separation per- 
pendicular to the satellite, p, was calculated from the physical distance between 
the receivers and the angle between the baseline and the incoming wave direction. 
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For cooparison with obaarvtstlons, a chaorctical analyala of the expcctad 
output aignal-to-nolaa ratio ia made in Appendix VI. Since the aatellite waa 
in motion, v* expected to obaerve a differential doppler aa a ramp in the phaae 
record. Appendix VII givea an example of the phaae change. 

3.2 RESULTS 

(a) The firat item of intereat la the differential phaae variance. Aa 
deacribed earlier, the phaae fluctuation von], be auperia^>oaed on a ramp because 
of the satellite motion. Samples 4^, 1 " l,...,n, are taken from selected 
portions of the chart recordings free from spurious disturbances. A least 
square linear fit was made for each run selected. Thus ^ ■ a -f k(t-t) 
vhere 



Results are presented in Table 3.1 for the runs made. It is seen that the runs 
phase difference fluctuation is less than one degree on the average. A typical 
chart record Is given in Figure 3. A and an average run is a»out 1 to 2 minutes. 

In computing the differential phase variance, we needed to consider the 
contribution due to random noise. It was assumed that it Is equal to 1/SKR, 
where SNR Is the slgnal-to-nolse ratio at the Input to the phase meter. 

Table 3.1 gives the observed SNR of the receiver at IF as estimated 
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rir’j-’.r 3.4 cjlvrt rfcorosr t;^\ce of Dirrr.RtNTiAi. p?l\sf. 

on the oscilloscope, rhe phase meter used, however, further narrows the 
effective bandwidth to 1 kHz from the 4 MHz in the receivers used. Thus 
there is an improvement of 36 dU over the receiver output. Towards the 
letter part of the table the SNTl was quite unequal for the two receivers 

because one of the receivers became noisier. The value given is the 
smaller of the two. The last column gives the standard deviation, after 

correction for the noise variance (“ • 

w>NK 

Assuming o? . “ • 

9, A f 9,N 

, is effective variance 

9»A 

is observed variance 
9 

c is the variance due to random noise 

it can be seen that the effect of noise is negligible for all cases since 

the effective SNR are all above 40 dB. 

(b) An attempt was also made to estimate power spectra from the dif- 

ftrc.-.tial phase records after removing the effect of the ramp. In the 
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ebsence of large sample sizes, the usual method of calculating the auto- 
correlation and Fourier transforming to the power spectrum cannot be used. 

The approach adopted is to compute the finite Fourier transform of the data 
samples using an FFT routine and then compute a smoothed periodogram, as sug- 
gested by Brilllnger [ 10] , where s(t) is an Integer with near \ and the 

averai^e is taken over (2m 1) adjacent perlodograms 



is the periodogram at frequency >. for a time series >:(t) of duration T. 

A theoretical basis of this technique is given in [ 10] 

In oar analysis, the value of m is chosen to be 2, i.e., the smoothed 
spectrum at each freque.icy A is on average of 5 adjacent perlodograms. 

The a’lalysis was made on several runs obtained; two typical spectra ob- 
tained are given in Figure 3.5. 

-5/2 

In both cases the power density spectra fell off as f . The contin- 
uous line through each plot is a rough average of the computed points. Also 
given is the 95 percent confidence Interval drawn as dashed lines. It is 
shown in [10] that 
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3.3 THEORETICAL INVESTIGATION 

A great deal of literature exlats on wave propagation through the 
Ionosphere and troposphere, t /eating them as random media with fluctuating 

dielectric constants. In the nonlonlzed lower atmosphere, the variation of the ^ 

dielectric constant Is due to winds, thermal currents, etc., while In the Iono- 
sphere, the variation Is due to deviations In electron density. Much work, e.g., 

[11,12,13], models the Ionospheric layer as a thin, phase screen. This approach 
allows considerable simplification In calculating statistical properties of the 
field. A second body of knowledge originates In Investigations of optical and 
millimeter wave propagation through the troposphere. Some review articles are 
[14,15,16], Experimental results have been reported in [17,18] for tropospheric 

propagation. Considerable effort has been given In applying the multiple- 
scatter technique to trans-ionospheric problems In the equatorial and auroral 
regions where strong scintillations have been observed [19]. It has also been 

reported, e.g.,[20l» that at radio freouencles In the GHz range the fluctu- 
ations in the v.ive are due primarily to the ionospheric irregularities. The 
most recent paper [22l provides both theoretical and experimental results of 
amplitude and phase scintillation caused by the ionosphere and appeared only 
after our experiment was completed. 

The objective of this study is to confirm the magnitude of phase dif- 
ference fluctuations observed in the experiment described earlier. 

At this latitude, the ionosphere is relatively quiet and since the satel- 
lite is in line-of-sight of the receivers, there is a strong specular com- 
ponent. Consequently the single scatter approach will be used, as treated 
[ 21 ]. 
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Following Tatarnkl (21], we shall assume the wavelength X Is much less 

than the inner scale of turbulence in order to simplify the calculations. In 

so doing, the scalar wave equation instead of the vector equation can be used. 

This assumption also means the effective scattering cone is small and results 

in considerable simplification. It has been shown in (23] that most of these 

results at optical frequencies can be extended to microwave frequencies. Start- 

_ llcz 

Ing with the wave equation, and F. ■ E^e $, for linear polarization, 

+ 7(E • Vfnc ) + k^e E ■ 0 , k^ ■ c (3.1) 

' r or o o o 

The -love assumptioi. reduces the problem to solving the scalar equation 


w!;er.. 

Also 




» 1 




f + 2ik-^ - 0 , k 


E “ E + E, , K, « E 
X ox lx lx ox 

A V +y, 

„ A . IS o 1 
E ■=■ Ae “ e 

X 


y " L, /E “ X + is 
1 lx ox 1 



(3.2) 


X and S are the log-auplltude and phase departure of the perturbation 
1 

term caused by scattering in the random medium. 

From (3*2)lt has been shown (section 7.3 of (21 J) that the phase 
structure function D^(a) can be derived i 


Ds(p)|j,“ E{lS^(x,y,7) - (x' ,y ’ , Z) | ^ } 

« An / U-J (kp)J F (x,o)xdK 
■' o o s 


(3.3) 

(3.4) 



where P 


is the antenna separation which is assuomd 
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to be less than L , the outer turbulent scale. In %rriting 3,4 we have 
o 

assumed locally homcgeneous fields. 


£(r) - f(0) + /// d-e^’' ’^)d«(ic) 


(3.5) 


where d^(tc) is a random amplitude satisfying 


E{ d^(»Cj^)d?* (< 2 )^ “ 
and also local isotropy . 


(3.6) 


D^(r) - E{]f(r+r^) - f(r^)|2} - D^(|rl) 


(3.6) 


Details are given in Tatarski's first book, chapter 1 and 7. In (3.4), the 
term F^(<,o) is the two-dimensional spectral density of the structure 
function of the phase fluctuation in the plane z » Z. As in Chapter 3 in 
the book, it is related to the refractive index spectral density 
f^(<^.K 2 , |z'-z"|, ^^— ) by 


Fg(<,o) “ ^ // cos [ ^ ^ J cos F^(K,|z'-z"|,^^)dz’dz 


2 / 7 _ _ » ' 


<3-7) 


is in general a function of the separation of the two points of interest, 
z'-z", and can vary along the path from o to Z. Thus assuming there are 
n layers each with a spectral density F^^, from (3.7) 


2 n-1 ziU 


F^OC.o) -X 


I II * cosi^), (3.8) 


i=o zi 


t - z'-z", n -(z'+z")/2 


z ■=• Z 
n 


0 , 
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This Is an extension of Che results in ( 21 J, Chapter 8, where Che irregu- 
larities are correlated only within each layer. Now ussuae can be 


written as 




o.») 


From (3.8)the double integral can be rewritten as 

. 2 . .2 


*i'*’*l+l . 


2(n-zj) 


,r 2 c^(r)dn / 


-2(r.-z^) 


i+1 , _ 


-2(n-z,„) 


"i'^*i + l 


2(n*ai^l) 


^ti^*'l’'^2’^'^ ^ ^ co8^^-jp=^) di 


As in Chapter S of [21] if <i i. 1, <^C/L << 1 and 2(h-r^), are 
of mat;nltude “ ’’t^' integrals of both terms in (3.10) are 
approy.inuitely 


2 1 1 + cos 


pJl)j / t^^(K-^,K-^,C)df. - 4n<^^(K) cos^[^-^p-> 


<: ^(<) is the 3-dinensional spectral density of the i medium so 

F (-:.o) - "V cos^[‘--g-^ dn 

i-o i 


(3.11) 


If hovever, C^(T) is constant within each layer, we can expand the inner 
n 

integral as in Chapter 7 of (21 J instead. 
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. 2 , ^2 


dC ■*■ *^°®V F^,(ie,o)dn 


Cl 


2 . .2 


/-I +Z ‘‘^ ^2 -C/2 * '°‘V f,/"'®) <•'' 


'i+1 i 


(z 




os^ (Z 2 ) 




(3.12) 


Equacion 3,8 reduces to 

F^«.o) . ^ I (1 + 


cos!^ (Z - ^ -|— ^ -) J 


(3.13) 


So bv (3.4) 


i^y(p) " J] ^ U*JqC<p)] ( ] »«dK (3.14) 

i»o o 


n-1 


where [ ] Is the tern (3,13) within square brackets. 

During its passage to earth, the wave passes through a layer of 
ionosphere, some free space and a thin troposphere about 10 km in height. 
As mentioned earlier, the tropospheric effect can be neglected unless the 
wave arrives at the receivers at a low’ elevation angle. Assuming we use 
(3.14) the structure function reduces to 
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D^(P) 


L / U-J^(kP)] (1+^ 8in!^ co*-'5^(Z-i) J «^(k)k<1k 

0 k^L 


(3.15) 


where <*^(*^) ■ 0 in the free space layer. 
We shall use the nodel given in (20 ] 

< c\ > 


’ 2 TTY 

o 


(3.16) 


where c^(r) - <c> (l+c^(r)), the relative permittivity. In the ionosphere 


- c (1-f ^/f^) - c (f 

Op o p N 


c e 
o r 


f 2 _ f 2 

. (r) „ _J1 — E) 

f o f 


<^T> 


.. <fM\2 > E 

n ^-'■'^2,2 
o ] 


(3.17) 


Ty;.ical values arc ’ 500m, = 7 Mhz, 


<(f)^> '■= 'O.l 

O 


Rather than evaluating (3.15) directly, we shall use Tatarskl's ap- 
proach hy calculating 0. (o) + D (p). It can he sho-vn, using (3.16), 

«\ s 


D,(o) + D,(p) - 2 tt k L / (1-J (<P)1 C (<)K'd< 

AS 0 O c 

“ 2<c 2> k^LL^ (1 -pK^(p/L )/L^l 


where is the modified Bessel function of second kind. 


(3.18) 
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m 

Da(p) - Fa(k,o)<«I< 

0 

- 2 <€?> kV [-^ eL V (p/L )-C(p) + C(o)l (3.19) 

^ AkL 4kL'* ° 

o o 

C(p) - / - J (*«P) sin ^ coi ^ 5TT 

o ° 2k L.k 2J 2 

o 

D (i3> is thus 3.19 subtracted from (3,18). 
s 

3.4 OBSERVATIONS 

Equation 3.19 is integrated numerically and a typical set of results 
is given in Figure 3.6 the parameters given there. The experimental re- 
sults all indicate very small phase fluctuations of less than one 

degree standard deviation even up to 800 ft separation. It implies very 
strong correlation of the signal phase. It is interesting that the expected 
gradual increase in phase fluctuation with baseline length was not observed. 
If we had an extra receiver so that the phase difference between a pair 
close together and the other far apart, simultaneous records would allow 
ua to make firmer conclusions. Alternatively if more time was allowed at 
the time of the experiment, the baseline could have been changed and the 
runb repeated for each baseline on the same day. 

The computed results using the model of the ionosphere given show a 
gradual increase in phase fluctuation with separation p. Note however that 
the electron density fluctuation of 10 percent was quoted in literature dis- 

cuv:.-r.g scintillations near the equator, which experiences much more activity 

than nt temperate regions. As given in 3.20, if <(^)2>** were reduced to 0.03 

N 



Phase std. deviation 
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froa 0.1 instead, the predicted fluctuations would be 10 times as small as the 
values In Figure (3.6) which are closer to the measured values. 

Tie frequency spectra of the phase difference appears to have slopes 

froc -2 to -2.5. In [ 16] > Tatarski derived the expressions for the power 

spectra but for the model we used as given in 3.15 and 3.16 the integral 

is not easily evaluated. But since the Kolumogorov spectrum used by Ta- 

tarshl goes as k while that in 3,16 is approximately proportional to 

-4 

K , the results are not expected to be too different. The theoretical 

—8/3 

power spectrum derived there predicts a frequency dependency of f which 
is quite close to that in Figure 3.5. A word of caution, however, is needed 
here. It must be remembered that the chart recordings are not strictly 
differential phase at an instant but contain a time factor owing to the plane 
wave arriving at an angle not perpendicular to the receiver baseline. The 
result for October 22 should be more realistic since the angle of arrival 
is about 5* with the baseline normal. The spectrum reported in ( 22] sh-ows 
a f ^ relationship. 

The differenticl phase fluctuations have been observed to be consis- 
tently small so that phase coherence can be guaranteed at least up to 800 ft 
in case of future experiments involving an extended random array. One possi- 
bility for future investigation is "o extend the baseline to several thousand 
feet perhaps using radio instead of wire links between receivers. i 

The investigation reported above Involves only the down-link situation 
where the wave incident on the atmosphere can be assumed to be plane. The up- 
link involves the propagation of a spherical wave into the atmosphere, which 
Fried has investigated for optical frequencies [24]. Except for different 
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valiMS, hi« fomulae should apply for tha Cranalonoapharic link Co a gaoaCaClon- 
ary aaCalllCa ualng microwave frequenclaa In analogy Co Che down-link case. 
Porccllo alao discussed sone aapecC of an uplink propagaCion [25]. 
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4. SUMMARY AND CONCLUSIONS 

It has been shown that the objectives of the project, namely the 
feasible design of an array in space and the carrying out of a spatial 
correlation experiment, have been fulfilled. 

For the system design part of the project, three radio camera system 
concepts have been described. The first one, based on the use of fixed 
references on the earth's surface to synchronize a randomly distributed 
antenna array which is in geostationary orbit, and then the array scans in 
the neighborhood of the synchronizing source, has been shown to be theoreti- 
cally feasible, but physically impractical . This is due to the requirement 
that the beacons can be uniformly located on any part of the earth, which is 
Impractical. 

The second system concept which was introduced and fully developed 
is based on the assumption that an element location subsystem can be created 
consisting of beacons clustered on the land masses and randomly distributed 
over the small angle subtended at the antenna array by the earth. An 
algorithm was introduced which showed this assumption to be completely 
valid. Thus phase measurements at each array element due to signals trans- 
mitted from the beacons are used to uniquely determine the positions of 
the elements of the array. Tolerance results were also given which showed 
the algorithm to satisfy all of the system specifications. Assuming the 
positions of the array elements are known, an algorithm which is closely 
related to the surveying algorithm is used to locate the position of trans- 
mitting sources. Further the electronic system also was designed, and It 
was calculated that this system concept is physically practical and it 
fulfills the desired objectives. 
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Finally a third ayaten concept, the Slaultaneoua Tranafora Array, was 

Introduces. No analysis was carried out; however It Is strongly felt that 

this new concept Is a possible alternative to the second system concept If 

subsequent design effort discloses problems not yet uncovered. 

The experimental part of the project was satisfactorily completed . The 

project plan called for a study of the spatial correlation of the phase of the 

received signal up to a distance of SCO feet. It was f.hown that even with up 

to 800 feet separation between the receivers, the standard deviation of the 

phase fluctuation was less than one degree. (This value Is so small that It 

Is regrettable that the project plan did not Include a much larger baseline.) 

The results were taken over a variety of atmospheric conditions and at different 

times of the day. A spectral analysis on the reduced data was also attempted 

and on those records which were computed, the phase power spectrum showed a 

-2.5 

frequency dependence of f ’ , wh^rh Is slightly lower than that published 

“3 

recently 122] which showed a f relationship. 
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5. RECOMMENDATIONS FOR FUTURE EFFORT 

The work to be done should be a combination of theoretical study • 
system design, and experimental verification of applications of large- 
scale, space-borne arrays. Specific topics to be considered are: 

(1) The mechanical design of large nonrlgld spaceborne array 
systems with capabilities of Internal station keeping. 

(2) Communication between the array modules and the space station 
(this Is dependent on the particular mechanical design chosen). 

(3) The design of up-llnk/down-llnk to Include radio connunicatlons 
from the earth to the ground source via the large array. 

(4) Antl-RFI adaptive nulling. 

(5) VFRC test of self survey and direction finding algorithms. 

(For this purpose It may be possible to study the results of the 
experiment being conducted In the summer of 1977 at VFRC In 
which the four elements of an array are surveyed using a set of 
beacons. The survey Is limited to the x-y plane only. It may 
be necessary to modify the experimental procedure to ascertain 
Its applicability to the space system design.) 

(6) Breadboarding and testing of an advanced array module for the 
spaceborne system. 

(7) Assembly of multi-element L-band, ground based array. Self 
survey from ground beacons. Use array to locate airborne targets. 

(8) Design the central processor. 

(9) Generation and distribution of electric power within the array. 

(10) Analytical study of the Simultaneous Transform Array and com- 
parative evaluation with the self-survey array system. 
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(11) Experimental design of ground (and sea) transmitting systems. 

(12) Experimental design of ground receiving system. 
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RELATIONS BETWEEN POINTING ANGLES FROM A SPACEBORNE ARRAY AND THE 
POSITIONS ON THE EARTH'S SURFACE [1.1] 

The steps involved in determining the longitude and the latitude 

of a point source radiating on the earth's surface from its direction 

cosines u and v measured in the array coordinate system are summarized 
a a 

in a concise form in Table I.l. 

In the first step we calculate the direction cosines of the radiat- 
ing source in the array auxiliary coordinate system fro® Its 

direction cosines measured in the array principal coordinate system 
(X Y Z ). Then, in the second step, we use the values calculated above 
to determine the direction cosines of the radiating source in the earth's 
auxiliary coordinate system (^^g^le^le^ * transformation of coordi- 

nates in the third step leads to the direction cosines of the radiating 
source in the earth's principal coordinate system (X^Y^Z^) . Finally, in 
the fourth step, one readily calculates the longitude and the latitude of 
the radiating jource. 

The array associated with the system of concern will be as large as 
40000 wavelengths. A design approach which assumes that the positions 
of the array elements (in the array coordinate system) are exactly known 
is too idealistic to be of practical value for arrays of such a size. 
Therefore we will allow uncertainties in the array element locations and 
will design the system accordingly, l.e., we will use reference sources 
located on the earth's surface to synchronize the array. Since determi- 
nation of the direction cosines apsumed by a radiating source in the array 

[I.l] T* A. Dzekov, "Relations between Pointing Angles from a Spaceoorne 

Array and Positions on Earth's Surface", VFRC QPR No. 18, August 1976, 

pp. 81-88. 
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X Y z : 

a a a 
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t, Y, Z, 
la la la 

Step 

II 

X Y Z 

X Y Z. 

la la^la 

le^le^le 

Step 

III 

X, Y, Z, 

X Y Z 

le le le 

e e e 

Step 

IV 








^ “ tan ^(— ) 

e u 

e 

* - 8ln"^(w ) 
e e 


TABLE I. 1-TRANSFORMATION OF COORDINATES FROM ARRAY’S COORDINATE 
SYSTEM TO THE EARTH’S COORDINATE SYSTEM 
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coordloace system requires knowledge of the direction cosines associated 
with the source which synchronizes the array, we will need a procedure by 
which the latter can be calculated. Such a procedure Is given in Table 1.2. 
The relations presented In the table are developed In a way analogous to 
that used to develop the relations In Table I.l. 

As seen, determination of the longitude and the latitude of a radi- 
ating source on the earth's surface requires hnowledg*! of the geocenter 
to array distance (H), the radius of the earth (R^) and the transformation 
coefficients* used In step I and step III In Table !•! or In step II and 
step IV In Table I. 2. The accuracy of position determination Is very much 
dependent on the accuracy by which the transformation coefficients are 
known. Precise knowledge of the transformation coefficients Is of sub- 
stantial Importance especially for a system which is expected to provide 
as high accuracy In position determination as that of the system we are 
concerned with. 

It Is unreasonable to assume that the orientation parameters of the 
array will coincide exactly with the parameters planned before launching; 
moreover it Is very likely that the orientation of the array will be time 
dependent. Hence, we must be concerned with the problem how one can 
calibrate the system, or In other words, how one can calculate (recalcu- 
late) the transformation coefficients when necessary. 

Calculation of the transformation coefficients Is equivalent to es- 
tablishment of the two auxiliary coordinate systems (X. Y, Z. and 

Xd XA 

^le^le^le^ the properties of which were described in the first part of the 

A 

The transformation coefficients In essence describe the orientation of 
the array with respect to the earth and the angular position of the ar- 
ray In the earth's coordinate system. 
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(H-R w, )/(R^ + - 2w, R H) 

e le e ie e 


V (1-w^ 
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u (1— 
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TABLE I.2-TRANSFORMi\TION OF COORDINATES FROM THE EARTH'S 

COORDINATE SYSTEM TO THE ARRAY'S COORDINATE SYSTEM. 
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discussion. Thsrsby, the accomplishments listed below and Illustrated in 
Figure I.3wlll suffice to obtain an adequate set of transformation coeffic- 
ients, l.e., to calibrate the system. 

(a) Determination of the direction cosines of the earth's center 

In the array coordinate system gives the orien- 

tation of the Z^^-axls in the array coordinate system. 

(b) Determination of the direction cosines of the array reference 

point in the earth's coordinate system. (X Y Z ); this gives the 

e e e 

orientation of the Z^^^-axis in the earth's coordinate system. 

(c) Determination of the direction cosines of an arbitrary point in 

space (the point can be located on the earth's surface) in both 

Che array's and the earth's coordinate systems. This point (the 

point A in the figure) together with the array reference point 

and the earth's center will define the plane Co be subtended by 

the Y-axis of the coordinate system and the X-axis 

of the X, Y. Z, coordinate system, 
le le le 



0 


E 




EARTH'S 

CENTER 


A (POINT ON THE 
EARTH'S SURFACE) 


FIGURE I.l UNIT VECTORS IN THE COORDINATE S YSTEMS OF 
THE ARRAY AND THE EARTH. 
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The treosformation coefficients will then follow from 


k, ■ s, a -f a, a -f a. a 
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k' ■ e, e + e, e + e, e 
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|k| X k’ 
k| X k^ 
|k| X k 2 


>2 — — — 
— - e, e + e, e + e, e 
lyx X lyy y lyz z 


— — X ki “ e, e + e, e + e, e 

1^, ^ I 1 Ixx X Ixy y Ixz z 


Therefore the only question Is how to obtain the components of the vectors 
k^, k^f kjl^ and In the corresponding coordinate systems. 

One way to determine the components of the vectors k^ and k^ In the 
array coordinate system is to Image the earth by a telescope mounted 
aboard the array. Namely, one can make use of the fact that known orien- 
tation of the telescope with respect to the array coordinate system and 
known characteristics of the optical system employed define a unique re- 
lationship between the direction cosines of the points In the object space 
in the array coordinate system and the position of the corresponding Im- 
ages in the image plane. Since the vector by definition points toward 
the earth's center its components can be calculated from the position of 
the earth's disc In the image plane. Concerning the vector k 2 we have In 
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prioclpl* ch« freedom co chooew ICe orientetlon (or the values of Ita 
componente in the errey coordinate aystem) arbitrarily. However, in order 
to facilitate the determination of the componente of the vector in Che 
earth'e coordinate ayetem (k^ by definition i.>aat lie in the aame plane aa 
k 2 ) the following procedure ahould be followed. We pick up an arbitrary 
object on the earth'a aurface (amall lake, email ialand, etc.) whoae lon- 
gltudu and latitude are known end whoae image can be identified in the 
L^age plane. Then, from the poaitlon of the object 'a image in the image 
plane we calculate the direction coainea of the object in the array co- 
ordinate ayatem. Finally, we require the vector k^ to be pointed towarda 
the object in queation, l.e., we consider the calculated direction cosines 
as cotponents of k^. 

Tl.e components of the vectors k| in the earth's coordinate system 
are specified by the longitude and the latitude of the array. (Methods 
exist by which the angular position of a geostationary satellite can be 
determined to within an rms uncertainty smaller than 0.1 railliradian [1.2]. 
This can be considered as quite satisfactory for the system we are con- 
cerned with.) Tl a components of the vector k^ are specified by the lon- 
gitude and the latitude of the object which was used to define the vector 
k^; this object plays the role of the point A in Figure 13 . 


fl 21 A.M. McCaskill, D.V. Neil, and A. A. Satterlee, "Launch and Orbital 
Injection of Intelsat IV Satellites", COMSAT Technical Review, 
Volume 2, Number 2, Fall 1972. 
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DEPENDENCE OF TdL UNCERTAINTY IN ANGULAR LOCATION ON THE BANDWIDTH 
OF THE SYSTEM 
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vill b« within ch« 3 dB region. It is easy to se« thst this region will 
have the font depicted by the shaded area in Figure II. 1. 



FIGURE II. 1-SYSTEM'S BANDWIDTH - SINGLE 
UNCERTAINTY RELATIONSHIP. 

Observe that the distance between the points C ari D corresponds to the 
array beanwidth and that the distance between the points A and B corres- 
ponds (and is aumerically equal) to the array bandwidth B^. 

What one can conclude from the figure is the following. The uncer- 
tainty involved in determining the angular coordinate of a source detached 
is an increasing function of the uncertainty within which the frequency 
coordinate of the source was knovm. Since in practice the uncertainty 
associated with the frequency coordinate is given by the system's band- 
width (3g* figure) one can also state that the uncertainty in 

angle is an increasing function of the system's bandwidth. Further, the 
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